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ABSTRACT We begin to investigate particle dynamics that is gov- 
^ ! erned by a nonstandard kinetic action of a special form. We are guided by a 

phenomenological scheme-the modified dynamics (MOND)-that imputes the 
mass discrepancy, observed in galactic systems, not to the presence of dark 
matter, but to a departure from Newtonian dynamics below a certain scale of 
accelerations, a a . The particle's equation of motion in a potential </> is derived 
from an action, S, of the form S oc Sfc[r(£), a a ] — f <fi dt. The limit a a — > cor- 
responds to Newtonian dynamics, and there the kinetic action must take 
the standard form. In the opposite limit, a a — > oo, we require — > 0-and 
more specifically, for circular orbits Sj~ oc a~ x -in order to attain the phe- 
nomenological success of MOND. If, to boot, is Galilei invariant it must be 
time-non-local; indeed, it is non-local in the strong sense that it cannot even 
be a limit of a sequence of local, higher-derivative theories, with increasing 
order. This is a blessing, as such theories need not suffer from the illnesses 
that are endemic to higher-derivative theories. We comment on the possibility 
that such a modified law of motion is an effective theory resulting from the 
elimination of degrees of freedom pertaining to the universe at large (the near 
equality a a ~ cH Q being a trace of that connection) . We derive a general virial 
relation for bounded trajectories. Exact solutions are obtained for circular or- 
bits, which pertain to rotation curves of disk galaxies: The orbital speed, 
v, and radius, r, are related by (the rotation curve) fj,(a/a )a = where 
a = f 2 /r, and fx(x) is simply related to the value of the kinetic action, S^, for 
a circular trajectory. Expressions for the angular momentum and energy of 
circular orbits are derived. Beyond the obvious matter-of-principle differences 
between the modified-law-of-motion theories treated here, and modifications 
of gravity, there are important applicative differences. For example, the ex- 
pression for the angular momentum-an invariant under temporal variations 
of a D -is different, and the evolution of disk galaxies, expected to ensue such 



variations, would be different in the two classes of theories. The various points 
we make are demonstrated on examples of theories of this class, but we cannot 
yet offer a specific theory that satisfies us on all accounts, even at the non- 
relativistic level that we consider here. We also explore, in passing, theories 
that depart from the conventional Newtonian dynamics for very low frequen- 
cies. Such theories may be written that are linear, and hence more amenable 
to solution and analysis than the MOND theories; but, they do not fare as 
well in explaining the observations. 

I. INTRODUCTION: THE GALACTIC MASS DISCREP- 
ANCY AND THE MODIFIED DYNAMICS 

It has been posited that the observed mass discrepancy, evinced by galac- 
tic systems, stems from the breakdown of the laws that are used to analyze 
their dynamics. This contrasts with the prevalent dark-matter doctrine that 
imputes the discrepancy to the presence of large quantities of, yet unobserved, 
matter. In particular, the problem has been attributed [1-3] to a departure 
from Newton laws that occurs in the limit of small accelerations, such as are 
found in galactic systems: There exists a certain acceleration parameter, a , 
much above which Newton laws apply with great accuracy. Much below a a 
the dynamics (of gravitational systems, at least) departs drastically from these 
laws, and may be epitomized by a relation of the form 

a 2 /a = MG/r 2 , (1) 

which gives the gravitational acceleration, a, of a test particle, at a distance r 
from a (point) mass M (instead of the usual a = MG/r 2 ). There are different 
theories that might spring from this seed relation. In fact, our purpose here 
is to point a new direction for searching such a theory. 

Arguably, the performance of this modified nonrelativistic dynamics 
(MOND) has been quite successful in reproducing the observations of galac- 
tic systems [4-7] (if not altogether without dispute-see e.g. Lake [8], and the 
rebuttal by Milgrom[9]). 

The following predictions are particularly clear-cut and well-nigh oblivi- 
ous to the details of the underlying theory. Most of them pertain to circular 
motion in axisymmetric disk galaxies. 

1. The speed on a circular orbit, supported gravitationally by a central body, 
becomes independent of the orbital radius for large radii [2]. This implies 
universally flat asymptotic rotation curves of isolated disk galaxies, in accord 
with high quality observations of extended rotation curves [4, 5, 7]. 

2. The asymptotic rotational speed, V^, depends only on the total mass, 
M, of the central body via = MGa a . This leads to a relation of the form 

oc L, where L is any luminosity parameter that is proportional to the total 
mass. This is to be compared with the observed strong — L correlation, 
known as the Tully-Fisher relation[10]. 

3. If one approximates elliptical galaxies by isothermal spheres, MOND pre- 
dicts [11] a strong correlation between the velocity dispersion in a galaxy, 
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a, and its total mass M: a 4 ~ MGa a , reproducing the observed Faber- 
Jackson[12] relation between a and the luminosity. 

4. MOND predicts[2] that galaxies with very low surface density-which 
correspond to low mean accelerations-will evince a large mass discrepancy. 
This has been forcibly born out by many recent observations of low-surface- 
brightness disks[13], and, with less confidence, for dwarf-spheroidal satellites 
of our galaxy [14]. 

5. According to MOND, elliptical galaxies-as modeled by isothermal spheres- 
cannot have a mean acceleration much exceeding a a (or a mean surface density 
much exceeding the critical surface density, E G = a G _1 ). This explains the 
observed cutoff in the distribution of mean surface brightnesses of elliptical 
galaxies, known as the Fish law[15-16] 

6. It is known that bare, thin disks, supported by rotation, are unstable 
to bar formation, and eventual breakup. MOND endows such disks with 
added stability that enables them to survive, provided they fall in the MOND 
regime, i.e. their mean acceleration is at a Q or below[17]. This can explain the 
observed marked paucity of galactic disks with mean surface density larger 
than the above E G , known as the Freeman law (see ref. [16]) for a discussion 
of this issue) . 

7. Above all, detailed analysis of rotation curves, for galaxies with high-quality 
data, show a very close match of the predictions of MOND with the observed 
rotation curves of these galaxies [4, 5, 7]. 

Generally, the dynamics of systems of galaxies-such as small groups, 
galaxy clusters, and the Virgo infall-which are not as well understood, can 
also be explained by MOND, with no need to invoke dark matter [3]. 

In the framework of MOND, all these tell us that a modification is called 
for when dealing with bound orbits. There is no observational information 
on unbound (scattering) trajectories, for massive particles. However, we may 
gain insight from observations of gravitational lensing of light from distant 
galaxies, by foreground clusters of galaxies[18]. These however concern rel- 
ativistic trajectories, which require a relativistic extension of MOND for a 
proper treatment (see below and Sec. IV). Henceforth we assume MOND for 
all types of nonrelativistic trajectories. 

The value of a Q has been determined, with concurrent results, by several 
independent methods; these rest on the twofold role of a a as the borderline 
acceleration (based on predictions 5 and 6 above), and as the scale of accel- 
eration in the deep MOND limit (based on predictions 2 and 3). The best 
determination comes from the detailed studies of rotation curves [4, 5, 7], and is 

a ~ (1 - 1.5)1(T 8 cm s~ 2 , (2) 

assuming the value of 75km/s/Mpc for the Hubble constant. Values of a/a a as 
small as 0.1 have been probed in studies of rotation curves; for the dynamics 
of the infall into the Virgo cluster one has typically a/a a ~ 0.01. 

MOND, which at the moment, is but a phenomeno logical theory, may be 
interpreted as either a modification of inertia, or as a modification of gravity 
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that leaves the law of motion intact [1]. True, General Relativity has taught 
us that gravity and inertia are not quite separable on the fundamental level. 
For instance, the Einstein field equations for the gravitational field imply 
the geodesic law of motion. However, at the level at which we work-i.e. a 
non-relativistic effective theory-it does make sense to distinguish between the 
interpretations. By the former we mean, foremost, a modification the must 
be implemented for whatever combination of forces determines the motion. 

In default of a full-fledged theory Milgrom[l,2,3] had initially used a work- 
ing extrapolating equation of motion of the form 

H(a/a )a = a N = -V0, (3) 

to calculate the acceleration a from the conventional potential <fi (®n is 
the Newtonian acceleration). The extrapolating function, /J,(x), satisfies 
//.(.x) x ~ > ' 00 ) 1 . so that Newtonian dynamics is attained in the limit of large 
accelerations. In the deep MOND limit we have fi(x) oc x, for x < 1 (and 
a a is normalized such that (i(x) ~ x), so that eq. (1) results, and all the 
above-mentioned predictions follow. Relation (3) is noncommittal on the 
matter of interpretation (see ref. [1]). As it is written it may be viewed as a 
modification of the law of motion for arbitrary potential <p, or, inverted to read 
a = — v(\ V0|/a o )V0, and assumed to apply only for gravitational potential, 
it may be viewed, instead, as a modification of gravity. 

Bekenstein and Milgrom[19] have then devised a nonrelativistic formu- 
lation of MOND that is expressly a modification of gravity: It replaces the 
Poisson equation (V • V0 = 4tcGq) by 

V • [M|V0|/a o )V0] = **Gq, (4) 

which determines the gravitational potential, <fi, generated by the mass-density 
distribution g. The acceleration of a test particle is then given by — V0. This 
formulation is derivable from an action principle, enjoys the usual conservation 
laws, and has been shown to be satisfactory on other important accounts[19]. 
It reduces to eq. (3) in cases of one-dimensional symmetry, and, in fact, even 
for less symmetric systems such as disk galaxies, has been shown to predict 
very similar rotation curves as eq. (3) . Notwithstanding the superior standing 
of eq. (4) , it is eq. (3) that has been used in all the rotation-curve tests of 
MOND[4,5,7]-being much wieldier. 

There are several relativistic theories with MOND as their non-relativistic 
limit [19] [20]; none is without problems. They are all based on modification 
of gravity, in that they can be formulated as a modification of the Einstein 
field equations for the metric, with particles moving on geodesies. Obviously, 
this entails an indirect effect on other forces when gravity is important; but, 
inasmuch as gravity can be neglected, they do not entail a modification of the 
motion of particles subjected to other forces. 

In this paper we begin to discuss formulations of MOND that, in contrast, 
may be viewd as a modifications of the law of motion of a particle in a given 
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potential field. The main implication here is that particle dynamics is affected 
for whichever combination of forces is at play. 

We derive the dynamics from an action principle, with an action S = 
Sp + Ski where the potential part, S p , takes the conventional form, while the 
kinetic action, Sj~, carries the modification. The latter is functional of the 
trajectory of the particle, and is constructed using the acceleration constant 
a a as the only extra parameter. The theory must approach the Newtonian 
dynamics in the limit a Q — > 0. When a a is not much smaller then all the 
quantities of the dimensions of acceleration, MOND takes its effect. The deep 
MOND limit corresponds to a a — > oo. To obtain the behavior epitomized by 
eq. (1) -or, to have asymptotic flat rotation curves for isolated galaxies-we 
expect Sk oc a^ 1 in this limit, at least for circular trajectories. 

Thus-in departure from the pristine statement of the MOND hypothesis- 
deviations from Newtonian dynamics appear, here, not only when the momen- 
tary acceleration on the trajectory is much smaller than a a . The special role 
of the acceleration is highlighted only through the introduction of the con- 
stant a a of these dimensions. On a circular, constant-speed trajectories, all 
the quantities with the dimensions of acceleration that can be built equal the 
acceleration itself. For such trajectories the acceleration only has relevancy 
for MOND (as is the case for rotation curves of disk galaxies). 

The main purposes of this paper are: (i) To adumbrate a possible ra- 
tionale for MOND, which we do in Sec. II. (ii) To demonstrate that it is 
possible to formulate MOND as a modification of the law of motion-in con- 
trast to a modification of gravity-and to expound some general properties of 
theories based on such a modification (Sec. Ill), (iii) To point out significant 
differences between the two interpretations, and to highlight the observational 
implications of these differences (Sec. VI). (iv) To present general, exact re- 
sults pertaining to circular motion, and point their relevance to the dynamics 
of disk galaxies. In particular, we give a general expression for the rotation 
curve of an axisymmetric disk galaxy, and derive simple expressions for the 
energy and angular momentum of circular orbits in Lagrangian theories; this 
we do in Sec. IV. 

Some of the points we want to explicate become more transparent in 
the context of a theory where the dynamics is modified in the limit of low 
frequencies, instead of the limit of small accelerations. Theories of this class 
may be constructed that are linear, and thus more amenable to analysis and 
solution; alas, as potential alternatives to dark matter in galactic systems, 
they seem to be at some odds with the observations. We discuss these, in 
some detail, in Sec. V, both for their own possible merit, and as heuristic 
auxiliaries. Finally, we list some of the desiderata that we have not been able 
to implement yet, and some questions that remain open (Sec. VI). 

We cannot pinpoint a specific theory. There is a large variety of theories 
that, as far as we have checked, are consistent with the available, relevant 
data (especially RCs). Also, we can offer no single model that possesses all 
the features we would ultimately require from a theory (see Sec. VI). 



- 5 - 



This approach opens a new course for searching for a sorely needed, rel- 
ativistic extension of MOND. 

II. PREAMBLE: THE POSSIBLE PROVENANCE OF THE 
MODIFIED DYNAMICS 

Various considerations lead us to suspect that MOND is an effective ap- 
proximation of a theory at a deeper stratum. Here, a most germane observa- 
tion[l] is that the deduced value of a Q [eq. (2) ] is of the same order as cH Q , 
where H is the present expansion rate of the universe (the Hubble constant) , 
and c the speed of light. (For H Q = 75 km/s/Mpc, cH Q ~7x 10~ 8 cm s~ 2 .) 
This near equality may betoken an effect of the universe as a whole on local 
dynamics, on systems that are small on the cosmological scale. 

There are, in fact, several quantities with the dimensions of an acceler- 
ation, that can be constructed from cosmological parameters [21], beside the 
above expansion parameter 

a ex = cH Q . (5) 

Another is 

a cu = c 2 /R c , (6) 
where R c is the radius of spatial curvature of the universe; yet another is 

a^c 2 \\\ l/ \ (7) 

with A the cosmological constant (A = A/3c 2 , where A is sometimes used). 

Only a ex is observationally determined with some accuracy, as the Hubble 
constant has been measured to lie between about 50 and 100 km/s/Mpc (see 
Huchra[22] for a review). On R c we have only a lower limit of the order of 
c/H 0: and on |A| we have only an upper limit [23-24] of order (H /c) 2 . So, 
only upper limits exist on a cu and a A , and these are of the order of a ex . In a 
Friedmannian universe, with vanishing cosmological constant, one has 

a cu /a ex = |1-0| 1/2 , (8) 

where O is the ratio of the mean density in the universe to the closure density. 
(We use these relations, derived from conventional cosmology, only heuristi- 
cally, as MOND may entail a modification of the cosmological equations.) So, 
if today is not very near 1, a cu and comparable. If the mean mass 

density falls much shorter of closure-as would be in keeping with the basic 
premise of MOND-but a non-zero cosmological constant renders the universe 
flat-as would conform with inflationary models [25] (see also ref. [24]) -we 
have today 

a x ~ a ex . (9) 

We see then that the near equality a a ~ a ex may be fortuitous, with a a 
being really a proxy for a c „, or for a a (or for another cosmic parameter). The 
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identification of the parameter behind a D would be an important step toward 
constructing an underlying theory for MOND; it is also of great impact in 
connection with formation of structure in the universe, and the subsequent 
evolution of galaxies, and galactic systems[21]: Whereas a x is a veritable con- 
stant, a ex , and a cu vary with cosmic time (in different ways). The matter has 
also antropic bearings. 

The Newtonian limit, which corresponds to a Q = 0, would be realized, 
according to the above identifications of a Q , when the universe is static (a ex = 
0), flat (a cu = 0), or characterized by a vanishing cosmological constant (ci\ = 
0), respectively. 

It is clear that to some degree the various cosmological parameters must 
appear in local dynamics. A finite curvature radius of the universe would en- 
ter, at least through the boundary conditions for the Einstein field equations, 
and will modify the local gravitational field. Likewise, a non-zero cosmological 
constant modifies the distance dependence of Newtonian gravity. The expan- 
sion of the universe also enters, for example by slowing down otherwise- free 
massive particles (or red-shifting free photons). All these are not, however, 
what we need for MOND. These effects are only appreciable over distances or 
times of cosmological scales-not relevant for non-relativistic galactic systems 
and phenomena. We thus do not expect H D , R c , and A to enter local dynamics 
as a time or length scale. On the other hand, a a is of the the same order as the 
galactic parameter of the same dimensions. So, if a a does remain as a vestige 
of the underlying mechanism, it will be the only such parameter important 
for non-relativistic galactic phenomena. In contrast, relativistic phenomena 
characterized by accelerations of order a Q or less, also involve scale length of 
cosmological magnitude. For example, a photon trajectory passing a distance 
r from a mass M, such that MG/r 2 < a a , has a curvature radius of cosmo- 
logical scale (see more in Sec. IV). So, the fact that only a a appears in MOND 
may be peculiar to the non-relativistic nature of this theory. 

The following analogy may help illuminate the situation we have in mind: 
Had all our knowledge come from experiments in a small, closed laboratory on 
the earth's surface, our dynamics would have involved a "universal" constant 
g-the free-fall acceleration. Unaware of the action of the earth, we would find 
a relation F = m(a — g) between the applied force F, and the acceleration a 
of a particle. Actually, g encapsules the effects of the mass, M e , and radius, 
i?0, of the earth. i? ffi , by itself, appears conspicuously in the dynamics only 
on scales comparable with itself; g alone enters the results of very- small- scale 
experiments characterized by accelerations not much larger then g. (The es- 
cape speed here plays the part of the speed of light in connection with MOND: 
g = vl s /2R & .) The effective departure here is of a gravitational origin, but we 
may describe it as a modification of inertia because it is in effect for whatever 
combination of forces F represents, even when there are no "local" sources of 
gravity (i.e. sources in addition to the earth). We would thus have no reason 
to consider it a modification of gravity. It may well be that the departure in- 
herent in MOND is also of a gravitational origin (and will then automatically 
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preserve the weak equivalence principle in the effective theory). 

But why, at all, should H Q , R c , or A enter local dynamics as an acceler- 
ation parameter? We can only offer a vague possible rationale: Even within 
the purview of conventional physics we expect a difference of sorts in the dy- 
namics of particles with accelerations above and below a , if a a stands for 
one of the above cosmological parameters. Consider, for example, a electric 
charge accelerated with acceleration a; the lower radius of its radiation zone 
is r rad = <? I a. To say that a <C a cu , for example, is to say that r rad is much 
larger than the radius of curvature of universe. In a closed universe, for exam- 
ple, this would mean that there in not even room for a radiation zone. Clearly, 
the radiation pattern must differ from the case where a ^> a cu for which the 
curvature is hardly felt within the near field. More generally, an accelerated 
particle carries a causal structure that includes an event horizon [26], whose 
scale is £ = c /a: This is the distance to the horizon; this is the space-time 
distance within which a freely falling frame can be erected; this is the typical 
wavelength of the Unruh black-body radiation, etc.. For accelerations that 
are low in the context of MOND (a <C a D ), the quality of the universe that is 
behind a a is felt within I. Conversly, we surmise that I is imprinted on the 
conjectured inertia field (e.g. by modifying the spectrum of the vacuum fields) 
in a way that is then sensitive to the acceleration of a test particle interacting 
with this field. 

All the above considerations may also usher MOND in as an effective 
theory of gravity. 

We thus envisage inertia as resulting from the interaction of the accel- 
erated body with some agent field, perhaps having to do with the vacuum 
fields, perhaps with an "inertia field" whose source is matter in the "rest of 
the universe" -in the spirit of Mach's principle. If this field bears the imprint 
of at least one of the cosmic parameters, an effective description of dynamics 
in which this inertia field is eliminated, may result in a scheme like MOND. 

We have not been able to put these vague ideas into a concrete use, and 
will have to make do, in the rest of this paper, with the study of effective 
theories. Nonetheless, we keep this conjectured origin before our eyes as a 
guide of sorts on what MOND may entail, beside the implications for the 
structure of galactic systems, which have been the phenomenological anchor 
for MOND. For example, such considerations would lead us reckon with the 
effects of an a Q that varies with cosmic time, as alluded to above. We may 
also suspect that if MOND, as we now formulate it, is only an effective theory 
delimited by cosmological factors, then it cannot be used to describe cosmology 
itself (as we would not describe the motion of earth satellites with a constant-fir 
dynamics in the above analogy). More generally, we noted above that MOND 
as described above-with a a as the only relevant parameter-may breakdown 
for any relativistic phenomenon with acceleration below a a (more in Sec. VI). 

There are attributes of the conventional laws of motion that we cannot 
take for granted in light of what we said above on the possible origin of the 
modified dynamics. For example, it is not obvious that the effective theory 
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describing MOND is derivable from an action principle. The question of which 
symmetries (conservation laws) the theory enjoys also becomes open. We shall 
certainly want to preserve the six cherished symmetries of space; to wit, the 
translations and rotations (or more generally the six isometries of a maximally 
symmetric 3-D space); what we said above raises no question with regard to 
these: the underpinnings of the effective theory are not expected to define 
a preferred position or direction. The same is true of parity. Regarding 
time reversal symmetry, if the expansion of the universe is in the heart of 
MOND, the fact that the universe has a well defined time arrow might well 
result in a strong TR asymmetry in the MOND regime, just as physics in an 
external magnetic field is TR non-symmetric. It would be very difficult to 
ascertain a breakdown of TR symmetry in the province of the galaxies. Many 
potential manifestations of this are forbidden by other symmetries, and there 
are other balking factors. The same is true of the time-translation symmetry. 
Here our guess would be that time translation might be affected through the 
introduction of the scale of time (H a l ), and this is not interesting, in the 
context of MOND, as it would lead to deviations only on very long time 
scales. 

Similarly, we are at sea regarding full Galilei invariance. While the ulti- 
mate underlying theory may well be Galilei (Poincare) invariant, or enjoy even 
a higher symmetry, this is not necessarily the case for the effective theory that 
results from the elimination of some degrees of freedom: the surmised "iner- 
tia field" , which comprises the hidden degrees of freedom, defines a preferred 
frame. Perhaps the effective theory does enjoy some generalization of Galilei 
invariance that involves parameter, and that reduces to Galilei invari- 

ance in the limit a a = 0. Imagine, for example, an accelerated observer in a 
curved space. As long as its acceleration is much larger than a cu of eq. (6) , so 
that its Unruh wavelength, or the size of the region affecting the dynamics, is 
much smaller than the radius of curvature, he will not detect departure from 
Lorentz invariance due to the curvature; for a <C a cu the non-zero curvature 
distorts the Unruh spectrum. 

As phenomenology does not yet require that we relinquish any of the 
above symmetries, and as we do not have any concrete theoretical impetus to 
do so, we shall retain them in the following discussion. 

We shall concentrate on effective theories of motion that are derived from 
an action principle, and conform to the MOND hypotheses. Such effective 
theories, which are obtained by eliminating some degrees of freedom from a 
more fundamental theory, are, in many cases, non-local; i.e., they are not 
derivable from an action that is an integral over a simple Lagrangian that is 
a function of only a finite number of time derivatives of the trajectory. A well 
known case in point is the Wheeler-Feynman formulation of electrodynam- 
ics[27], which results from the elimination of the electromagnetic field, and is 
non-local when only the particle coordinates remain as degrees of freedom. We 
shall show-without bringing to bear the origin of MOND- that an acceptable 
theory of MOND must, in fact, be strongly non-local. 



9 



III. GENERAL MOND ACTIONS 

One starts by assuming that the motion of a non-relativistic test particle, 
in a static potential field 0, is governed by an equation of motion of the form 

ml=-V0, (10) 

where m is the mass of the particle, and A (replacing the acceleration a in 
Newtonian dynamics) depends only on the trajectory r(t) (not on the po- 
tential); its value is, in general, a functional of the whole trajectory, and a 
function of the momentary position on the trajectory, parameterized by the 
time t. Of A we further require the following: (i) That it be independent of 
m, so as to retain the weak equivalence principle. In what follows we take a 
unit mass for the particle, (ii) According to the MOND assumption, A may 
be constructed with the aid of the acceleration constant a a , as the only dimen- 
sional constant. This is assumed to be the case in the non-relativistic limit, 
(iii) A has to approach the Newtonian expression, a, in the limit a a — > 0, for 
all trajectories, to achieve correspondence with Newtonian dynamics, (iv) The 
deep MOND limit corresponds to a a — > oo; phenomenology points to A oc a~ 1 , 
in this limit. We do not know that this shoul hold universally (i.e. for all sort 
of trajectories etc.), we shall assume this behavior for circular orbits, and for 

— * 

general trajectories assume only A — > 0. We also do not have observational 
information on the very limit a a — > oo; reliable information exists only down 
to a/a Q ~ 0.1. Here we shall assume this everywhere for large a a , but most of 
what we say is independent of this assumption, (v) We assume a a to be time 
independent, and so we are free to assume that A does not depend explicitly 
on time. If a Q actually varies with cosmic time, our results may be considered 
the lowest order in an adiabatic approximation. Explicit time dependence 
could also enter A if the particle has variable mass, which is simple to include 
(but which we bar), (vi) A is invariant under translations, and transforms as 
a vector under rotations of the trajectory, (vii) We do insist on full Galilei 
invariance. 

Indication that we need require some invariance connected with a veloc- 
ity boost could have come from the study of unbound trajectories. When 
a fast, unbound particle of velocity v scatters gravitationally off some mass, 
with pericenter distance r, the characteristic acceleration around pericenter is 
a ~ fg S /r, where v es is the escape speed from r. We can, however, with the 
aid of v, construct other quantities with the dimensions of acceleration, e.g. 
(in = (v n 4jpn(i) 1 / ( - n+1 \ Around pericenter, a n ~ (v /v es ) 2n ^ n+1 ^a^ which may 
be much larger than a itself. This might lead, in theories such as we discuss 
here, to much less of a mass discrepancy for such unbound particles, compared 
with that manifested by bound particles; for the latter, all quantities of the 
dimensions of acceleration are of order a. It is the appearance of the velocity v 
that engenders such terms. Galilei invariance would bar the appearance of v in 
the equation of motion, and all the permissible quantities with dimensions of 
acceleration are of the order of the acceleration, near pericenter. We may then 
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expect that bound, and scattering test particles will point to mass discrepan- 
cies of the same order (a /a). Observations of lensing by galaxy clusters may 
illuminate this question (for these v/v es is typically a few hundred), pending 
the understanding of the relativistic extension of MOND. 

While, for convenience, we treat motion in a static potential, all our 
results apply to the relative motion in the center-of-mass frame, of a system 
of two particles , interacting through a potential that depends on the relative 
position. So Galilei invariance of A for the particles assures invariance for the 
system. 

Equation (3) satisfies all the above requirements, but does not conform 
to another the we now add; viz., that the equation of motion be derivable 
from an action principle. As far as we can see, there is nothing amiss, in 
principle, with effective equations of motion that are not so derivable. If their 
symmetries do not lead to conservation laws, this may just reflect interaction 
with the hidden degrees of freedom. In fact, exactly such effective theories 
have been proposed as models for dissipative systems[28]. That being as it 
may, our purpose here is to investigate equations of motion with an action 
underpinning. This is a strong constraint that enables us to say much about 
the theory. One result of this assumption is that a theory whose kinetic action 
is Galilei invariant, and has the correct Newtonian limit, and the required 
MOND limit, cannot be local; i.e. A in eq. (10) cannot be a function of a 
finite number of time derivatives of r(t). (From this we understand why eq. 
(3) , which is local, cannot be derived from an action.) Other results concern 
a virial relation, conservation laws, and adiabatic invariants. In yet another 
we shall show that circular orbits are necessarily solutions of the equation 
of motion in a potential with the appropriate symmetry, and derive explicit 
expressions for the rotational speed. 

We formulate the action principle as follows: There exists a kinetic action 
functional Sf c [f(t),t 1 ,t 2 ,a ] (not necessarily Lagrangian) that is a functional 
of the trajectory, r(t), given between two end times t ± , and t 2 - Further assume 
that the variation of S). under a change Sr(t) may be written as 

5S k = —T~ l f * A(t) ■ 8r(t) dt + T" 1 ^ • 8r{t)% (11) 

where A does not depend on the end points. The second term on the right- 
hand side is (up to the T normalization) the difference between the values of 
S • 5f(t), at the two end times; S is some linear vector operator, acting on 
functions of t, [and may depend on r(t)}. We then consider only equations of 
motion (10) whose kinetic part A can be obtained in this way from an action 
functional. The normalization by the time interval T = t 2 — t 1 is introduced 
for convenience, so as to render the action finite when t — > oo; it does not 
affect the equation of motion. 

All the above requirements we made of A are now carried, mutatis mu- 
tandis, to the kinetic action S^. 
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Because we want the kinetic part of the equation of motion to bear the 
modification that MOND entails, we choose to discuss the kinetic part sepa- 
rately. The full action of the theory is 

S = S k + S p , (12) 

where the potential part has the usual form (except for our normalization): 

S p = —T~ l / </> dt. (13) 

The variation of the full action is 

SS = — T _1 f 2 [A(t) + V0] • 8f(t) dt + T" 1 ^ • 8f(t)]ll. (14) 
J tl 

Note that we do not define the physical trajectories as extrema of the 
action; that would be the case only if one restricts himself to variation, 8r(t), 
for which the ends term vanishes, which is usually done by fiat. The trajec- 
tory may well be an extremum, but we need not and wish not make such 
requirements. 

Standard Lagrangian theories are special instances (see below), but the 
above definition will suit us in more general cases. When the theory is non- 
local, it usually makes sense only when we take ti to — oo, and t 2 to +oo. In this 
case the equation of motion is a non-local equation for the full world-line of the 
particle. A finite trajectory (specified between finite times) is permitted if it 
is a segment of a world-line that solves the equation of motion. An important 
question regarding a given theory is how many (vector) parameters are needed 
to specify a permitted world line; this is also the number of initial conditions 
that one has to specify for a finite time t. 

Consider now variations of S under coordinate transformations. For so- 
lutions of the equation of motion, only the ends term remains in eq. (14) , 
and it is thus clear that S generates conserved quantities, when underlying 
symmetries of S are present. For example, if S is invariant under translations 
f — > r + 5e, we see that the difference of the values of (9 acting on 1) 
between the two end points must vanish. This quantity may be identified as 
the momentum. Similarly, if S is invariant under rotations about a direction 
n, since Sr = en x r, we find that — n ■ [9f x f(t)] is the same at the two end 

— * 

points. Thus we identify — S x f(t) with the angular momentum. 

Now take a dilatational increment in a solution of the equation of motion: 
8f(t) = ef(t). For trajectories on which S • f(t) does not increase as fast as t, 
at t — > ±oo we get, from eq. (11) , for the variation of (for the full world 
line T — > oo) 

8S k = -e(A-r p ), (15) 
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where () is the time average over the trajectory. This certainly holds for 
periodic trajectories, but surmisedly for any bounded trajectory. We now 
calculate SS k differently: On dimensional grounds 

S k [\f(t),^o] =X 2 S k [f(t),a ], (16) 



for a constant A. Putting A = 1 + e, and identifying the first order terms in e, 
we get, as in Euler's theorem for homogeneous functions, 

SS k + ea ^ = 2eS k . (17) 
oa a 

Comparing with expression (15) for 5S k we get 

-(A.r) = 2S k -a ^, (18) 
and from the equation of motion 

5 fc (l + ^) = ^(r-V0), (19) 
where Sk = -§£-4^. This is a generalization of the virial relation (in the 

CJLiL do 

standard Newtonian case S k is the mean kinetic energy, with our normaliza- 
tion). 

More generally, if S k depends on a number of dimensional constants, 
C e , £ = 1, i", with dimensions (length)"? x (time) 13 ?, respectively, then S k 
in eq. (19) is replaced by — J2l=i a e d\n ZT e ■ we re pl ace ^fc with 5", and put 
cj) = we get a more general relation that reduces to the former when S = 
S k + 4> with <j) depending only on r, because then we must have <j) = C^f(f/r ) 
with C and r of dimensions of velocity and length, respectively.) 

In Newtonian dynamics A = a is linear in r(t). In a MOND theory this 
may not be: From dimensional arguments 

A[\r(t),\a ] = \A[r(t),a ], (20) 



or, taking A = a a 1 , 

A[f(t),a ]=a A[a; 1 f(t),l}. (21) 

If A is linear in f(t), it is independent of a a altogether. 

It may be more convenient, sometimes, to write the equation of motion in 
terms of the expansion coefficients of f(t) in some basis of complete functions, 
instead of r(t) (which is itself the expansion coefficient of the trajectory in the 
functions 5(t' — t)). This is not a common practice because (i) for time-local, 
Lagrangian theories the locality is manifest only when t is used as an indepen- 
dent variable (ii) t is a natural choice when the potential is time independent. 
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If we forgo locality anyway, the benefits of using a different decomposition of 
the trajectory may outweigh the drawbacks. For example, we may use the 
Fourier transform, p(ui), of f(t). As a functional of p the kinetic action may 
take a particularly simple form. We discuss examples below. 

It is well known that theories with higher derivatives tend to suffer from 
various diseases, we discuss these and their prevention at the end of this 
section, after we have considered a few examples. 

LAGRANGIAN THEORIES 

We now specialize to kinetic actions that are derivable from a Lagrangian 
function: 

S k = T~ l L k dt (22) 

(the total Lagrangian being L = L k — (p). We assume <p to depend only on r 
(and, in particular, not to depend explicitly on the time). In this case A is 
given by the standard Euler-Lagrange expression 



N 



-1 



' dL L 



l(<) 



(23) 



and the operator S is 



N 



(24) 



with 



N 



Vr, 



9L k 



- (i—m) 



(25) 



Here, D = j^, and a superscript (i) stands for the ith time derivative. 

If the number, N, of derivatives on which L k depends is finite, the theory 
is local. Even when it is non-local, the theory may be a limit of a sequence of 
local ones with Lagrangians , N — > oo. (A common example is an Lj, that 
can be expanded in a power series in some parameter-a Q in our case-where the 
coefficients are function of a finite number, N, of derivatives.) When it is not 
even that we call the theory strongly non-local. For such theories the sums 
in the expressions above do not converge for at least some trajectories. These 
expressions may still be of value, as they may converge for some trajectories, 
and may be analytically continued for even more trajectories. A different 
expression for A must then be given, to establish an equation of motion. 

As explained above, translation invariance implies the conservation of the 
momentum 3(1), which here is just p*i (this can also be seen from the fact that 
A = ^jjr)- (In a many-body system with bodies of masses m i5 and potentials 
that depend on coordinate differences, J2i m iAi = 0.) Similarly, the angular 



n 



momentum-whose component along n is conserved when is symmetric under 
rotations about n-is given by 



N 



N 



N 



m=l 



or 



JV 



m=l 



9L k 



x 



£=1 m=l 



(i— to) 



- (i—m) 

i 

(26) 
(27) 



This also follows straightforwardly from the fact that r x A = so the 
rotational symmetry of <fi [which reads V0 • (n x r) =0] implies, from the 

equation of motion, d ^/^ = 0. 

As in standard dynamics, if the potential is rotationally symmetric about 
n, n-J is conserved, even when the Lagrangian depends explicitly on time, such 
as when a Q varies with cosmological time. Obviously, this has implications for 
galaxy evolution that is driven by such variations. 

If L does not depend explicitly on time, the energy H is conserved, where 



N 



-l^i-m+l) . 



i = l m = l 



dL 



(m-l) 



L. 



(28) 



We separate H = H k + <fi, and call H k the kinetic energy. The conservation of 
H can be deduced by noting that v ■ A = H k + so that, from the equation 
of motion, if = 0, we have H k + <fi = 0. 

— * — * 

If we add to a total derivative to L, then A, H, ,p ± , and J remain the 

— * 

same, but p { for i > 1, and hence 3 do change. 

In appendix B we derive an expression for the energy in general La- 
grangian theories, which, specialized to the case of MOND, reads 



H = L k -2a ^ + <P + Q, 
da n 



where 



N i-1 



g = -^^(-l)- 1 (z-m)f 



=<i-m) 



9L k 



(m-l) 



(29) 



(30) 



This has an immediate corollary: for any solution of the equation of motion 
on which Q does not grow as fast as £, at t — > ±oo, such as on periodic, and 
presumably other bounded trajectories, we can write 



H = (H) = (L k ) - 2a Q 



da r , 



+ (0>- 



(31) 
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Combined with the virial relation (19) this gives 

(# fe > = ^7<r-V0>, (32) 

where 7 = (1 + 2S k )/(l + S k /2). In the Newtonian limit 7 = 1; in general 
7 < 4; for circular orbits in the MOND limit 7 = 2. 

If the theory is Galilei invariant and is local or weakly non-local, L k may 
be written (up to a total time derivative) in the form 

L k = l -av 2 +L k (a ^ 2 \^\...), (33) 

where a is a number, and L k depends only on the acceleration and higher 
derivatives of f(t). We show this in appendix A. Correspondence with Newto- 
nian dynamics requires that a = 1, and that L k vanish for a a = (or reduce 
to a time derivative, which we can eliminate from the outset). On the other 
hand, the MOND limit requires a = 0. Thus, Candidate, Galilei-invariant 
theories for MOND, that are derivable from an action, must be strongly non- 
local. This is the basis for our statement that eq. (3) , which is local, and 
Galilei invariant, cannot be derived from an action. 

Higher-derivative Lagrangian theories that are not strongly non-local may 
be cast in a Hamiltonian form as described by Ostrogradski[29]. Very suc- 
cinctly, one defines N coordinates, the first of which is gl = f; the rest are 
defined recursively by q m = m = 2,...,N. There are N momenta, 

which are those defined by eq. (25) . The Hamiltonian is H, as defined in 
(28) , in which the derivatives of f(t) have been expressed by the 2N coordi- 
nates and momenta, taken to be independent. The Hamilton equations of the 
usual form then give the equation of motion, and, to boot, the above relations 
between the coordinates. The Hamiltonian formalism is particularly useful 
when one endeavors to quantize the theory; we shall have little use for it here. 

SOME HEURISTIC MODEL ACTIONS 

Before embarking on general issues concerning properties of the solutions, 
we describe a few examples of theories on which we shall demonstrate various 
points. 

Consider first a class of kinetic Lagrangians of the form 

L k = \vf{a%-'D)v. (34) 

<^= =>■ 

Here, D and D are the time derivative operators acting to their left and right 
respectively, and g is a function of the derivatives of f(t) from the second up- 
ward, having the right dimensions (l 2 /t 6 ) to make the argument of / dimen- 
sionless. We must have /(0) = 1 for Newtonian correspondence; furthermore, 
as stated in point (iii) in the beginning of this section, for L k to approach the 
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Newtonian limit for all trajectories, f(z) has to be non-singular at z = in 
the complex plane. If it has a Taylor expansion, with a non-zero radius of 
convergence: 

oo 

f(z) = l + J2 b mZ m , (35) 

m=l 

then 

1 ' 1 ' d ^\o(„, % ,36) 



1 1 1 

T 2,- L i. 2-12,- L i. 4 

fc = 2^ 2 ^ a 2 ° 



eft 



which has the form (33) . Clearly, may not be represented by this ex- 
pansion for all values of a a , lest it not be strongly non-local, and not have a 
MOND limit. In the next section we derive an exact expression for the rota- 
tion curve in this class of theories, and show that to obtain asymptotically flat 
rotation curves we must have f(x) oc x~ 1/2 for x — > oo (x real). For example 
f(z) = (1 + z)~ l/2 satisfies all the above requirements; it has a finite radius of 
convergence. 

The Lagrangian (34) is Galilei invariant: the required Newtonian limit 
implies that f{a 2 Dg^ 1 iy) acting on a constant vector v , to its left or right, 
gives v itself; so, under a transformation f(t) — > r(t) + v t + r*o, changes 
by Vq/2 + v ■ v(t), which is a total time derivative. 

Second, consider Lagrangians of the form 

L k = \vf[<{DhD)-']v. (37) 

Here, Galilei invariance is achieved because of the required MOND limit 
/(oc) = 0. 

We consider also actions of the form (22) , where depends also on 
parameters, 5^, that are functionals of the whole trajectory (and not functions 
of the local values of r(t), and its derivatives). Take, for instance of the 
form (34) [or (37) ] with g (or h) of the form 

g={G)=T- 1 J Gdt, (38) 

where G is a function of the derivatives of f(t). The action is not, strictly 
speaking, Lagrangian; we keep that epithet but call such Lagrangians com- 
pound. When all the g^s are of the form g i = (Gj), with local G;S, the equation 
of motion has the same form as the Euler-Lagrange equation, with replaced 
by 

M k = L k + J2Gi^. (39) 

In an even wieldier class of (non-Lagrangian) theories we define the kinetic 
action % as a solution of the equation 

s* = ^-' [g.F (40) 
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The integrand is defined when v(t) can be written as a sum of exponentials 

exp(at), with arbitrary complex a: e at F (^ Sk ^ D ^j e pt = F (%^) e (a+/3) *. In 

the Newtonian limit F(x) x ^°° > 1. To get the correct MOND limit for circular 
orbits, we must have F(x) = (x/t]) 1/3 for x <C 1 (with our normalization of a Q , 
j] = 9/8). F(z) is analytic at infinity, according to our general requirement, 
but is non-analytic at z = 0. This theory is "almost" linear, whence come its 
heuristic value. 

Note that only trajectories that can be expanded in real frequencies are 
relevant here: For an infinite-time world line there is no finite-^ solution to 
eq. (40) , if there are complex- frequency components to f(t). Such orbits 
necessarily approach the Newtonian regime, for large (positive or negative) 
times, but they are irrelevant there as candidate solutions of the equation of 
motion. 

The equation of motion is 

(1 + S k /2)F (-^-) D 2 f(t) = -V0 (41) 



with S k = 

To obtain eq. (41) we write 



/ * v(t)F' (^pj f At). (42) 

The two terms in the first integral are equal, as DD acts symmetrically. The 
second integral can be written in terms of S^., leading to 

SS k = (1 + SJ2) 1 - fdt^-F m . (43) 

Now, when f(t) and 5f(t) can be expanded in Fourier components with real 

frequencies, the last integral equals — f dt 5r(t) ■ F ^— ^ k J? j D 2 r(t), up to 

end terms [when F(z) is analytic at z = this is seen by integration by parts; 
when we go via Fourier space the end terms are lost] , from which follows eq. 
(41). 

To insure Galilei invariance, L k must actually be invariant, not just in- 
cremented by a time derivative; S k itself must be invariant-not just changed 
by end terms-as can be seen from the equation of motion. This, in turn, is 
insured by the required MOND limit [F(0) = 0]. (Hence the choice of DD 
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in the argument of F, in the definition of S^, is essential. A choice of — D 2 
instead would have given an expression that differs by an ends term, and that 
is not exactly invariant under Galilei transformations.) We want to be 
determined uniquely by eq. (40) . If this has many solutions, for a given 
trajectory, we need an added prescription to choose from among them. The 
obvious desideratum that Sf~ go to its Newtonian value, for that trajectory, 
when a a — > may do (see next subsection); it eliminates Sj~ = 0, which is 
always a solution. Note that the theory is time-reversal symmetric (if is a 
solution so is its value for the time reversed orbit, and the two have the same 
Newtonian limit). 

Finally, we bring for comparison a local Lagrangian theory with 

Lk = \v 2 f{al/a>); (44) 
obviously it not Galilei invariant. 

THE NATURE OF THE SOLUTIONS OF NONLOCAL THEORIES 

Higher-derivative theories are notoriously problematic. Pais and Uhlen- 
beck[30], and more recently, Jaen et al.[31], Eliezer and Woodard[32], and 
Simon[33] review these problems in Lagrangian theories, describe possible 
solutions, and give further references. Very succinctly, the problems are: (i) 
The energy as defined in eq. (28) is not bound from below: there are solu- 
tions of the equation of motion with arbitrarily small (large negative) energy, 
(ii) There appear runaway solutions of the equation of motion, on which the 
velocity, acceleration, etc. blow up at t — > oo. In linear theories the diver- 
gence is exponential, (iii) When N is the highest derivative on which the 
Lagrangian depends, the equation of motion is of order 2N ', and require the 
specification of that many initial conditions. While this may not be a matter- 
of-principle problem, it is disconcerting. There is a connection between these 
three problems, and they can usually be eliminated in one fell swoop, by 
adding appropriate constraints to a finite order theory, or by culling out the 
majority of the solutions of the equation of motion that are not analytic in 
some expansion parameter [3 1,32, 33]. 

We can see the relation between the runaway character of a solution and 
the energy from eq. (31) : When the solution is not runaway, the assumption 
underlying this relation is valid, and then we see that H is bound from below, 
since Lf~, and— are positive definite. 

When the theory is of infinite order, we may expect exacerbation of the 
above problems. For instance, one may have to dictate an infinite number of 
initial conditions, and this may signify a total loss of the initial-value formu- 
lation [32]. 

Can MOND-oriented theories avoid these daunting problems? As has 
been emphasized by Eliezer and Woodard[32] a strongly non-local theory (in 
our terminology) is not a higher-derivative theory, and need not suffer from 
any of the problems endemic to such theories. We saw that Galilei-invariant 



19 



MOND theories must be strongly non-local and may well be trouble free. 
Theories that emerge as effective approximation to healthy theories, by an 
elimination of some degrees of freedom, are obviously as healthy as their pro- 
genitors. As explained in Sec. II we envisage MOND to come about in just 
this way. 

We now demonstrate that the theory whose action is defined by eq. (40) 
, with a (solvable) harmonic-oscillator potential (p = ^uj 2 r 2 , is free of the 
above problems. As explained below eq. (40) , only trajectories that can 
be expanded in Fourier components with real frequencies have a chance of 
solving the equation of motion (41) . From this equation we also see that 
every frequency u> that appears in the expansion must satisfy 

(1 + S k /2)F u 2 = < (45) 

We can choose F to give a unique (positive) solution for cu 2 , given Sj~. The 
physical trajectories have thus a single frequency, and can be written as 

r -(t) = -L( r >^ + r -e-^). (46) 

The defining equation of S^, (40) , which supplements the equation of motion, 
reads here 

Sk = \^r%F , (47) 

where = r • r|J. We can also write from eqs. (45) (47) 

{l + S k /2)S k = l -rlul (48) 

Take for instance 

F{x) = {l + r ] x- 1 )- 1 '\ (49) 

which has the desired Newtonian and MOND limits (with 77 = 9/8). It is easy 
to ascertain that i?e[f ] and Im[f ] determine u 2 (and S k ) uniquely. 

Thus, there are no runaway solutions, and the physical trajectories are 
determined by two initial conditions, just as for a Newtonian harmonic oscil- 
lator. Here too the motion is harmonic, but the frequency depends on the 
amplitude. For example, in the deep MOND case-i.e., when u 2 r <C a G , we 
have S k rs 1, and F(x) {x/rf) 1 ^, and uj is determined by 

-1/4 

> u 2 . (50) 

A circular orbit is a special case which we shall discuss below. 

The case of an anisotropic harmonic potential, which is as easy to solve, 
is also instructive. The motion is still harmonic with a unique frequency in 
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each of the principal directions of the potential, the equation of motion now- 
being 

(l + S k /2) F (^-yi=u 2 0m , m= 1,2,3. (51) 

with 

3 

(l + ^/2)^ = ^|r 0m |^ 2 m - (52) 

m=l 

Given the three amplitudes and phases encoded in the complex amplitudes 
r 0m one determines and the three frequencies uj m ; each depends not only 
on its own amplitude; they are coupled through S^. 

In general, exponentially exploding solutions are precluded in MOND, 
because they posses derivatives that all blow up at infinite time, and this 
necessarily takes us to the Newtonian regime where they are not physical 
trajectories. In finite-order MOND theories even trajectories that runaway as 
a high enough power of t are precluded. This, in itself, does not solve any 
of the problems listed above. We saw that the criterion for a non-runaway 
solution is that Q of eq. (30) not grow as fast as t at ±oo. Take, for instance, 
the local MOND Lagrangian given by eq. (44) . Clearly, the equation of 
motion cannot have solutions for which a increases indefinitely at t — > oo; 
still, the theory has all the usual illnesses: The runaway solutions are such 
that v oc t 1 / 3 for t — > oo, so a — > 0, but the Q term in the expression for H is 
oc v 3 , enough to qualify as a runaway. 



IV. CIRCULAR ORBITS: DYNAMICS OF AXISYMMETRIC 
DISK GALAXIES 

No other observation is as central to the discussion of the mass discrep- 
ancy in galactic systems, and as clear-cut, as the measurement and analysis 
of rotation curves of disk galaxies. Fortunately, the relevant trajectories are 
particularly amenable to treatment in the framework that we consider here. 
These are circular, constant-speed orbits (the latter epithet will be implicit 
hereafter) in the mid-plane of an axisymmetric, and plane symmetric potential 
field, 4>. 

THE EQUATION OF MOTION-ROTATION CURVES 

We now show that, under the MOND assumption, and for a general action 
of the type given by eqs. (11) - (13) , the speed, v, on a circular orbit, in the 
mid-plane of a potential with the symmetry of a disk, depends on the orbital 
radius, r, through a relation of the form 

^{a/a )a = (53) 

where fx(x) depends on the specific kinetic action, and a = v 2 /r. Thus, eq. 
(3) gives the exact expression for the rotation curves, and existing rotation- 
curve analyses, which have all employed this expression[2,4,5,7], are valid in 
the present framework. 
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For consistency we first ascertain that, for a circular orbit, the functional 
derivative, A, in the equation of motion, is a radial vector of a constant mag- 
nitude. This is a property of the kinetic action alone (irrespective of the 
potential) and is necessary if circular orbits are to be solutions of the equation 
of motion in the mid-plane. This property follows from the assumed trans- 
lation and rotation invariances of the kinetic action. Under a infinitesimal 
translation r — > r + e, S k changes by e ■ f A dt, according to eq. (11) for 
the variation of S k [the ends term vanishes as f(t) and df(t) are periodic]. 
This must vanish for arbitrary e, so J A dt = 0. From the symmetry of the 
orbit, the component of A perpendicular to the orbit must be constant and 
hence must vanish. Similarly, the invariance of S k under a rotation about an 
arbitrary axis n (5f ccnxr] implies that J Ax f(t) dt = 0. Thus, A cannot 
have a tangential component either, and is a radial vector of constant mag- 
nitude, as required. The action underpinnings is essential; it is easy to write 
an equation of motion (10) in which A has the above symmetries but is not 
radial for circular orbits. (But TR is then not a symmetry.) 

— * 

The magnitude of A depends only on the parameters of the orbit, to wit 
v, and r. Because A/ a is dimensionless, it must, by the MOND assumption, 
be of the form fj,(a/a ), leading to 

A = /u(a/a )a, (54) 

and to eq. (53) , in turn. For there to be a circular-orbit solution for every 
positive value of xfx(x) must take all positive values, and A must be di- 
rected inward. Uniqueness requires that xji{x) is monotonic. As a corollary 
we get that any kinetic action that is rotation and translation invariant nec- 
essarily gives the Newtonian relation for circular orbits, if it does not contain 
any dimensional constant. 

There exists a simple and useful expression for fj,(x) in terms of the values 
of the action for circular orbits: 

»(a/a ) = 2v- 2 S c k (l + ±S c k ), (55) 

din S c 

where Sf, is the action calculated for a circular orbit, and Sf, = — ^ a k . This 
relation enables us to circumvent the calculation of functional derivatives. 

Equation (55) follows from the general expression (18) for (A ■ r)-which 
holds for any periodic orbit-with eq. (54) , and the equality a ■ r = —v 2 . 

With our normalization, S k has the dimensions of velocity square; so, 
must be of the form 

S c k = l -v 2 \(a/a ). (56) 

(Any functional of the orbit reduces to a function of v, r, and a , for a circular 
orbit. Thus, S?/v 2 is a dimensionless function of these, and must be of the 
form \\{v 2 /ra a ).) Hence, 

/*(*) = \(x) + ^ = X(x)[l + X(x)/2], (57) 
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where X(x) = din X/dln x. 

Only the time- reversal-symmetric part of S k contributes to S%, and hence 
to //, because, for a circular orbit, time reversal is equivalent to some 7r-rotation 
about an axis in the orbital plane, under which is invariant. 

When the action is derivable from a Lagrangian , L^-simple or compound, 
local or not-that does not depend on the time explicitly, one has S% = L^, 
where L| is the constant value of the Lagrangian on the circular orbit. Then 

v(a/a ) = 2v- 2 L c k (l + h c k ), (58) 

.,, f c din LI 

Consider, for instance, Lagrangian actions with L k of the form (34) . The 
function g is constant on a circular orbit, and, in fact, must be proportional to 
r~ 4 i> 6 , on dimensional grounds. The proportionality constant may be absorbed 
into a a . The operator DD, when sandwiched between any two derivatives of 
r c (t) (or — D 2 , when acting on any of these), is just the multiplication by 

J2 = v 2j r 2_ We then find that 

L% = \v 2 f{al/a% (59) 
from which results the expression 

V(a/a ) = f(x)[l + f(x)]\ x=a 2 /a 2, (60) 

where f(x) = —din f /din x. 

Since /i does not depend on the choice of g, many theories give the same 
rotation curve for a given mass distribution. For similar reasons, Lagrangians 
of the form (37) give the same rotation curve (with the same choice of /), 
although they are very different theories. The same is true of the compound 
Lagrangians discussed below eq. (37) , and of the local Lagrangian given by 
eq. (44) . 

For actions of the form (40) , we find from the general expression for 
/i(x) eq. (55) [or directly from the equation of motion (41) ] that /i(x) = [1 + 
t}\(x)]\(x) , where X(x) is the solution of X(x) = F[x 2 X(x)/2] with the correct 
Newtonian limit; i.e., the solution, for a = 0, of X a (x) = F[x 2 aX a (x)/2] that 
satisfies X a (x) — > 1 for a — > oo. 

Clearly, we can add to the Lagrangian a time-reversal-antisymmetric part 
without affecting the rotation curve. The same is true of a more general type 
of terms that vanish on circular orbits [e.g. oc (v ■ a) 2 ]. 

The fact that the rotation curve, for a given mass distribution depends 
only on the value of the action for circular orbits is gratifying in that it allows 
us to test this general class of theories by rotation-curve analysis, without 
having to know the exact theory [by adopting a best value for a a , and a 
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reasonable form for fJ,(x)]. This, however, leaves us bereft of an important 
observational guide in discriminating between theories within the class. 



THE ENERGY AND ANGULAR MOMENTUM 

When the action is derivable from a Lagrangian, and energy and angular 
momentum are defined-in the sense that the formal sums that express them 
converge, or may be understood as analytic continuations-we find simple ex- 
pressions for these quantities on circular orbits. 

Because scalar functions are constant on circular orbits, Q in expression 
(29) for the energy vanishes and we have 

H c = L c k -2a -^ + <j>. (61) 

In terms of X(x) of eq. (56) , which now describes the dependence of Ljr 
on a/a a , we can write 

H c = ^v 2 u(a/a ) + <j>, (62) 

where 

v(x) = \(x)[1 + 2\(x)]=»(x)k(x), «(g) = 1 + 2A(g) (63) 

1 + X(x)/2 

In the deep MOND limit, where fx(x) ~ x, and X(x) ~ 1, u(x) ~ 2x. 

Putting together eqs. (62) , (63) , and eq. (53) [fi(a/a )a = <j>'] which 
governs circular orbits in a disk potential cj) (taken to vanish at infinity), we 
can write for such a motion 

H c = -n(a/a )v 2 (K + 2/4>), (64) 

where 4> = din 4>/dln r. In Newtonian dynamics (// = 1, k = 1) there are no 
negative-energy circular orbits for potentials with <p < —2. It is interesting to 
note that in the strict MOND limit we have k — > 2, and the limiting potential 
approaches the Newtonian or Coulomb potential, with <f> — > — 1. 

Coming now to the angular momentum, we note that, for a circular orbit, 
vectors in the orbital plane [f(t) and its derivatives, derivative of with 
respect to those, etc.] are radial or tangential, according to whether they are 
even or odd under time reversal, respectively- because, for such orbits, time 
reversal is equivalent to some it rotation about a diameter of the orbit. We 
see then from expression (27) for J that only the time-reversal-symmetric 
part of L/j contributes. Furthermore, the summand in eq. (27) can now be 
written (after successively shifting time derivatives from the second factor to 
the first) as r (m_1) x ■ Hence, we can write for a circular orbit 

N f)T 

m=l 
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Since r i - m l) = —uj 2 Q x f( m) , with Q = r 2 rxv,we have 



N 



J c = u~ 2 u y mr» • (66) 



Finally, from expression (B-4) in appendix B, for this sum, we get 



2 / dL c \ 

J e= - 2 ^L t -a^y (67) 

In terms of X(x) of eq. (56) we have 

J c = V (a/a )f x v, v (x) = X(x) [1 + \(x)} = ^(x) 1 + Hx) (68) 

1 + A(x)/2 

So, in the very MOND limit J c — > 4w 3 /3a , and, interestingly, depends only 
on v. 

In modifications of gravity, the expression for the angular momentum is 
the same as in conventional dynamics; so, the adiabatic invariance of J has 
rather different consequences in the two classes of theories based, respectively, 
on modifying gravity, and modifying the law of motion. 



V. MODIFIED DYNAMICS AT LOW FREQUENCIES 

Notwithstanding the phenomenological success of the MOND scheme, it 
is of interest to consider modifications that hinge on other system properties. 
After all, galactic systems differ from laboratory and planetary systems in 
attributes other than the typical accelerations in them. Here we examine 
modifications that hinge on the time scale involved. More accurately, we have 
in mind theories that involve a constant, O, with the dimensions of frequency. 
In analogy with the above discussion of MOND, the equation of motion is of 
the form 

A[f(t),Q] = -V0. (69) 
From dimensional considerations we find that 

A[f(t),n] = 2 A[f(Ot),l], (70) 

and that 

A[Xf(t),Q} = XA[r(t),Q], (71) 

Thus, A is necessarily homogeneous in r(t), and there is nothing to prevent 
us from having A that is linear in r (t) . Theories with a linear kinetic term are 
rather more amenable to analysis; in particular their equation of motion is easy 
to derive from the action, when they are thus derivable; hence their heuristic 
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value. Henceforth we shall limit ourselves to such equations of motion, whose 
general form is 

f(D/n)D 2 f(t) = -V0. (72) 

In the present context the appropriate Newtonian limit is O — > 0, if we 
want approximate Newtonian behavior for laboratory phenomena, where the 
typical frequencies would be much higher than O. Thus, 

(73) 

Galilei invariance of A is assured if f(x) is not singular at x = 0. 

If we wish to follow this course beyond heuristics, and propound such 
theories as potential alternatives to dark matter, we have to note the following: 
For circular orbits in the mid-plane of a disk potential to be solutions of the 
equation of motion f(x) must be an even function of x: f(x) = h(—x 2 ). Then, 
eq. (72) gives for such orbits 

h[(v/rn) 2 ]v 2 /r = ^. (74) 
Flat rotation curves for disk galaxies result if 

%) ~ V 1/2 , for y « 1. (75) 

The asymptotic rotational speed, V^,, is then related to the total galactic mass 
by = MGVt. This is at some odds with the observed Infra-red Tully-Fisher 
relation (V^, oc L IR ), if the IR luminosity of a galaxy is proportional to its 
mass (see point 2 in Sec. I). The value of O would have to be about !0~ 3 H o to 
account for the mass discrepancy in galaxies (this value is approximately the 
orbital frequency of the sun in the Galaxy). Typical dynamical time scales 
in clusters of galaxies are much larger than in galaxies, so we would predict 
rather higher discrepancies in clusters. Not much of an Oort discrepancy is 
expected in the solar neighborhood, as the frequencies of the perpendicular 
motion are at least a few times higher than Q, and due to the linearity of the 
kinetic term there is no direct effect of the low-frequency orbital motion on 
the perpendicular motion (there is only a shift in the equilibrium state for the 
oscillations). This linearity also assures the desired non-Newtonian center-of- 
mass motion of composite systems (stars, clouds, or atoms in galaxies, etc.) 
even though their components rattle inside them at high frequencies. 

Phenomena in the laboratory or solar system are typified by frequencies, 
uj, much higher than O, and if the correction to Newtonian laws, in the limit 
uj O, is of order (Q/c<;) 2 -as expected if f(x) is even in its argument-it would 
not be detectable with present experimental capacity. 

The evenness of / in its argument is tantamount to time-reversal (TR) 
symmetry of the equation of motion, and also to its derivability from a La- 
grangian action with 

Lk = \v{t)h[^]v{t)- (76) 
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(Only the time-symmetric part of the Lagrangian affects the equation of mo- 
tion, which is thus necessarily TR symmetric.) 

Pais and Uhlenbeck[30] have discussed in detail the properties of a class 
of theories with Lagrangians equivalent to that of eq. (76) , with h(z) an 
entire (integral) function. They show that such theories suffer from the prob- 
lems discussed in Sec. III. Our insistence on the particular Newtonian limit 
(73) (which they are not concerned with) puts us outside the province of their 
analysis: An entire function has an essential singularity at infinity (barring 
the conventional case of a constant /), and thus cannot have the Newtonian 
limit as required in eq. (73) . (For a valid Newtonian limit, f(z) must go to 1, 
as z — > oo, in all directions in the complex plane, if the Lagrangian is to ap- 
proach the Newtonian Lagrangian for all trajectories. For example, harmonic 
trajectories correspond to an imaginary argument of /; time-exponential tra- 
jectories correspond to real ones, etc..) Our theory must then be non-local, 
indeed strongly non-local. In the present context this follows only from the 
required Newtonian limit; in the case of MOND, Galilei invariance, and the 
MOND limit were also brought to bear on the analogous conclusion. As we 
explain in Sec. Ill, the fact that our theory must not be entire in the time- 
derivative operator is a blessing: only such theories have a chance of being 
free of the difficulties that characterize higher-derivative theories. 

VI. DISCUSSION 

Examining modified laws of motion in the framework of the MOND 
scheme, we find that if the underlying kinetic action has the correct Newtonian 
limit, has the phenomenologically inspired MOND limit, and is Galilei invari- 
ant, the theory must be strongly non-local. It is this class of theories that we 
have concentrated on. They are not infinite-order limits of higher-derivative 
theories, which are known to have severe problems. We have demonstrated 
that MOND non-local theories may be formulated that are free of these prob- 
lems. An attractive possibility is that MOND results as a nonrelativistic, 
small-scale expression of a fundamental theory by which inertia is a vestige 
of the interaction of a body with "the rest of the Universe" , in the spirit of 
Mach's principle, or with an agent such as vacuum fields. The approximate 
equality of the acceleration "constant" of MOND, a G , and cH Q , may be a 
witness. Such effective theories are usually strongly non-local. 

We are not able yet to pinpoint a specific effective theory. The study 
of rotation curves-which provides the most detailed and clear-cut evidence 
for a mass discrepancy-is a strong test that is passed by all theories of the 
genre we consider here; but, it can only constrain the values of the kinetic 
action on circular orbits. Further constraints could come from the detection 
of mass discrepancies in light bending. Beyond these, existing observations 
of the mass discrepancy in galactic systems are hardly helpful in discerning a 
superior candidate theory. Additional discriminating power might come from 
the study of systems with non-circular motions, such as elliptical galaxies and 
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clusters of galaxies, because theories that predict the same circular motion 
usually disagree on non-circular trajectories. 

Also potentially helpful is the study of a possible mass discrepancy man- 
ifested by the motion perpendicular to the plane in disk galaxies, which is a 
small perturbation on the circular motion. There may be such a mass dis- 
crepancy observed in the solar neighborhood, in the Milky Way (the Oort 
discrepancy), but the issue is still moot. Different theories of the present class 
predict different degrees of such a discrepancy. For example, it was pointed 
out [2] that eq. (3) predicts a discrepancy of magnitude 1/ fi(a c /a ) where a c is 
the acceleration in the circular orbit (as long as the acceleration perpendicular 
to the disk is small compared with a c ). Analysis of small-perturbation, per- 
pendicular motion in the theory defined by the action (40) , gives a different 
value for the discrepancy: + S?/2)F(S?oj 2 /a 2 )], where u is the frequency 
of the perpendicular motion. Had uj been the frequency of the circular motion 
(actually <C uj), this would have given l/ / u(a c /a Q ). 

The examples we discuss are still formulated as test-particle dynamics. 
We cannot yet offer an example of a theory that fully accounts for the correct 
center-of-mass motion of composite objects: A composite body (e.g. a star), 
with constituents whose internal motion is Newtonian should still undergo a 
center-of-mass motion similar to that of a test particle, which may be well in 
the MOND regime. This we deem only a technical obstacle; a feeling that is 
supported by the following observation. We may formulate a MOND action 
as an integral over the Fourier transform, r w , of the trajectory, such that, 
in an extreme case, different frequencies are decoupled from each other, and 
for each frequency the motion depends on the value of uj 2 r UJ /a Q . Thus, the 
internal motions for which this value is large, obey approximate Newtonian 
dynamics, while the center-of-mass motion, characterized by a small value, is 
test-particle like. The theories discussed in Sec. V have this desired property. 
For various reasons we do not endorse such total spectral decoupling of the 
motion in the kinetic action (as occurs in the Newtonian action). 

The action in eq. (40) does couple different frequencies, and also give 
a partial solution to the problem: Consider a trajectory of a constituent of 
a "star" that performs "galactic" motion with (say one) frequency uj g and 
frequencies uii that describe the internal motion. We see from eq. (40) that, 
if the internal velocities are small enough compared with the galactic one, 
the value of is determined by the "galactic" component, even if uji uj g 
because F saturates at a value of 1 for large arguments. The condition for 
that is 

( Vi /v g y « S k uj 2 g /a 2 . (77) 

Then, the "galactic" motion will be test-particle like, and the internal motions 
will depend on the value of (Sj c /v 2 )(a 2 al), where a; = v^. This theory does 
not provide a fully satisfactory solution because it does not cover the case 
with internal velocities that are not smaller than the external one (motions 
within the atoms or nuclei are also internal). It does indicate, however that 
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the solution is within reach. 

Another desideratum we have yet to implement is an extension of the the- 
ory to the relativistic regime: a matter-of-principle must, that is also needed 
to treat light bending. We feel that understanding the underpinning of MOND 
is a prior goal, and it achievement will surely bring about the understanding 
of relativistic MOND. We note in this connection that for a phenomenon to be 
both relativistic, and to involve accelerations of order a a or below, it must also 
involve a length scale of cosmological magnitude. For instance, for the acceler- 
ation near the Schwarzschild radius of a mass M to be at a Q or below we must 
have this radius no smaller than the Hubble radius. Similarly, light trajecto- 
ries in regions where the Newtonian acceleration is of order a a , has a radius 
of curvature of the order of the Hubble radius. Thus, the relativistic-MOND 
regime falls in the large-scales regime. In light of the surmised cosmological 
origin of MOND this means that a relativistic extension of MOND might not 
simply hinge on the single acceleration constant a Q , which may, in fact, loose 
its role altogether. 

Like all effective theories, ours must have validity bounds that stem from 
the approximations underlying the theory. These delimiters are usually not 
recognized from within the theory itself; but, sometimes, inconsistencies or 
other problems are vestiges of these approximations. Consider, for instance 
trajectories with constant velocity; in all the model actions we have studied 
the correspond to the extreme MOND limit in which the action vanishes. 
Clearly they do not have a Newtonian limit, which is one of the axioms of our 
treatment; we see a hint of a delimiter here. We should also not be dismayed 
if our effective (relativistic) theory fails to account for cosmology; as explained 
in Sec. I we certainly expect a delimiter there. 

Finally, we list some differences between modifying gravity, and modifying 
inertia, beyond the fact that the latter pertains to arbitrary combinations of 
forces. We saw that conserved quantities and adiabatic invariants are different 
in the two approaches. In gravitational MOND theories, any isolated mass 
produces a confining gravitational potential: To obtain asymptotically flat 
rotation curves we must have an asymptotically logarithmic potential. It is 
not clear to us whether there a general answer to the analogue question in our 
class of theories (i.e., are there unbound trajectories about an isolated mass?). 
The approximation underlying our theory may, in fact, break down for such 
unbound orbit, as indicated above. 

APPENDIX A 

We derive the general form of a Galilei-invariant, local, kinetic La- 
grangian, [which can be written as a function of f(t) and a finite number 
of its derivatives] . We make use of the commutation rules 



and 



(A -2) 



[These stem from Q = (v-j^ + a- -^j + ...)Q.] For the equation of motion 
to be invariant under translations it is necessary (and sufficient) that V r £fc 
is some time derivative Taking the r-curl of this equality we get, with 
the aid of (A-f), that ^[V r x Q] = 0. Now, when a quantity F depends 
on a finite number of derivatives of f(t), and jjjF = (i.e. F is constant 
along any trajectory), then F must be an absolute constant, independent 
of any of the variables: If F depends on r 1 - * 1 , r^, r ( - k ' > only, and ~j^F = 
• ^) + ... + f (fe+1) • ®_^ k) = 0, for any choice of r*°\ f* k+1 \ then, clearly, 
all the partial derivatives of F must vanish. (Seen differently by noting that 
any two finite sets of derivatives are connected by some trajectory, so F must 
take the same value on any two sets.) Thus, V r x Q is a constant of the theory, 
and hence must actually vanish, because of the isotropy of the theory, which 
forbids the construction of a no n- vanishing, constant vector. Thus, Q = V r q 
for some g, and hence V r (£fc — q) = 0, and is independent of r up to the 
time derivative <j, which we can discard. 

Take then Lj, that depends only on the derivatives of r(t). For the equa- 

tion of motion to be Galilei invariant we must have VyL^ = for some Q. 
The left-hand side does not depend on the variable r, so take the derivative 
with respect to r to obtain ^[V r <S> Q] =0. As above, the tensor in paren- 
theses must be a constant, and is dimensionless, so it must be aSij, with a a 
dimensionless constant. Thus, Q = ar + U, where U is independent of r. We 
then have V v L>k = otv + ^jjr. Taking the v-cnrl (which now commutes with the 

— » — * — * 

t-derivative on U) we find that V v x U is a constant, and hence must vanish, 

— * 

and hence U is a f-divergence of some u that does not depend on r. From all 
this we get 

V v (L k -^av 2 -u) =0. (A - 3) 

The expression in parentheses is, thus, some that does not depend on for 
v as stated in eq. (33) . 

If Lfc is not of finite order, but is a limit, for iV — > oo, of a sequence of 
Lagrangians, , of finite order N, then, each of the 's is of the form (33) 
, and, hence, so is (if the LjJ 's are not themselves Galilei invariant, we 
can take their Galilei-symmetrized parts, and they too go to for iV — > oo). 



APPENDIX B 

Here we derive expression (29) for the energy function, which applies 
whenever this function can be defined. Starting with the general expression 
for the energy, eq. (28) , and shifting time derivatives from the second factor 
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to the first, we can write (for m > 2) 

dL 



^i-m+l) _ 



0fW 



(m-1) qt 

= q l + (-l) m - 1 ^ 



with 



m— 1 



= -E(-d' 



1 = 1 



dL 



<9f« 



(B-l) 



(S-2) 



Upon taking the double sum in eq (28) , the first term in eq.(B-l) gives 
Q with Q given by eq. (30) (after some rearrangement and summation). 
The second term in eq.(B-l) gives £\ zf^ ■ ^§y, which we now evaluate. 
Imagine that L depends on a number of dimensional constants Cg, £ = 1, I, 
with dimensions (length) ae x (time) 13 *, respectively. The Lagrangian itself has 
dimensions (length) 2 x (time)~ 2 . 

Considering the redefinition of the units of time we have 

L[X^C U X^Cj, A- 1 ^, X-^\ ...] = X~ 2 L[C U Cj, r^, r< l) , ...]. 

(5-3) 

Taking the derivative with respect to A at A = 1, we find, as in the Euler 
theorem for homogeneous functions, 



dL 



N 



IT 



=<0 



-2L. 



(5-4) 



e=i c i=i 

Putting eq.(B-4) together with the above, we get 

dL 



H = L + J2faC e — + Q 



i=i 



(5-5) 



We separate the Lagrangian into kinetic and potential parts L = L k — <fi; <f>, 
which is taken to depend only on r, must be of the form Cf(r/r ), where C has 
dimensions of velocity squared, and r those of length. Thus the contribution 
of (p to the sum in (B-5) is 20, and none to Q. We can then also write 



dLi 



H = L k + J2^C e ^ + <p + Q, 



i=i 



(5-6) 



and specializing to the case of a single constant a a ((3 = —2), we get eq. (29) . 

Another useful relation is obtained by considering, instead, redefinition 
of the units of length, which gives 

L k [X^C u A tt 'C„ ArW, Ar« ...] = X 2 L k [C u C u r« ...], 

(5-7) 

from which follows, as above, 



JV 



<=1 



dC, 



i=l 



dL k _ 

dm ~ k - 



(B-8) 
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